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ABSTRACT 

The extraction of the photodissociation cross sections of from 
Goulomb dissociation experiments is investigated. A carefnl stndy is 
done on the contribntions of the El, E2 and Ml mnltipolarities to the 
breaknp. A comparison with the data of a recent experiment is per¬ 
formed. It is shown that the extraction of the radiative captnre cross 
sections ^Re(p, 'j)^B which are relevant for the solar nentrino problem is 
not affected appreciably by Gonlomb reacceleration. A non-perturbative 
model is used for the purpose. Emphasis is put on the perspectives for 
future experiments which are planned at the University of Notre Dame, 
RIKEN (Japan), and GSI (Germany). An analysis of the total yields of 
“photon-point” processes in inelastic electron scattering is also done. 
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1 Introduction 


The solar neutrino problem is related to a factor of two difference between the the¬ 
oretically calculated and the experimentally measured neutrino flux from the Sun. 
This discrepancy is more accentuated for the measurements of the Homestake and 
the Kamiokande |0 experiments. The Gallex and SAGE ^ experiments also im¬ 
ply a discrepancy, although more moderate, between theory and experiment. While 
the Gallex and the SAGE experiments measure neutrinos which are as low in energy 
as 233 keV, the Homestake and the Kamiokande experiments measure high energy 
neutrinos. Thus, the discrepancy seems to be larger for the higher neutrino energies. 

Most of the high neutrinos {E,^ > 2 MeV) come from the beta-decay of The 
main mechanism through which the is produced in the Sun is the '^Be{p, 'j)^B 
reaction [|I|. If the cross section for this reaction were half of the value which is 
presently used as an input to the Standard Solar Model the Homestake and the 
Kamiokande experiments would agree well, although the gallium results would re¬ 
main unexplained. It is also worthwhile to say that this change might have other 
implications which can be in disagreement with the SSM model. A proposed solution 
to the solar neutrino problem is that the discrepancy lies not in the nuclear cross 
sections that are involved in the determination of the neutrino fluxes, but rather in 
the behaviour of the neutrinos during their flight to the terrestrial detectors, most 
likely matter induced oscillations into other neutrino types (The MSW effect [|^) to 
which the detectors are insensitive. However, not all the nuclear cross sections are 
known with the desired accuracy, and their inaccuracies will significantly affect the 
interpretation of the Davis ^'^Cl (Homestake) and Kamiokande Gherenkov detectors. 
The reaction '^Be{p, 'j)^B accounts for the bulk of the neutrinos seen by ^'^Cl and is 
the sole source of neutrinos seen by Kamiokande. 

The '^Be{p, 'y)^B reaction occurs in the Sun at a temperature corresponding to 
a relative p-^He energy of about 20 keV. Due to the Goulomb barrier, the magni¬ 
tude of the cross section at this energy is too small to be measured directly. Direct 
measurements have been done down to an energy of 120 keV [0. The measurements 
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of this cross section have always been difficult, because the target is radioactive and 
relatively short lived (jBe has a halflife of 53 days), and the measurements have been 
plagued with normalization. Thus, although the most recent and careful experiments 
of Fillipone et ah 0 have a quoted error of 10%, Barker and Spear argue that 
the true result may be smaller than that used in the predictions of Bahcall [^, by as 
much as 30%, due to a combination of normalization, data selection and excitation 
function uncertainties. 

Another controversial point is that the input in the solar models is an extrapolation 
of the cross section to the 20 keV region. This procedure yields an astrophysical S- 
factor close to 20 keV.b |^. These extrapolations are done based on nuclear models 
which allow some flexibility, and only the experiment can tell how accurate these 
models are. Therefore, as long as a measurement of the radiative capture cross 
section for the reaction ’^Be{p, 'y)^B at low energies is not done, the solar neutrino 
problem cannot be accepted as solved. 


A new technique []^ is now available for the measurements of these cross sections. 
One studies the inverse reaction ®i? + 7 — Be + p through the Coulomb breakup of 
on a heavy target. The Coulomb held of the heavy nucleus provides the equivalent 
photons. The Be{p, j^B cross section is obtained by detailed balance. This method 
has been shown to work very well for at least two reactions |^, 0 (for a recent review 
see |T^). Experiments with real photons are infeasible. But, in hrst-order. Coulomb 
excitation involves the same matrix elements as those of photo-induced processes. 
Furthermore, the coherent held of a large Z nucleus results in very large dissociation 
cross sections, which can be studied in present experimental facilities |^. A recent 
experiment by Motobayashi et ah have used the Coulomb dissociation method 
to extract the photodissociation cross section ^B + j — Z Be -|- p. 

The main problems with the method are (a) to experimentally separate the nu¬ 
clear from the Coulomb contribution in the breakup events, and (b) contribution 
from higher order ehects, usually called by reacceleration, or post-acceleration, ef- 
When the nuclear contribution to the breakup cross section is 


fects [|1^ [^, 

large compared to the Coulomb breakup mechanism for the experimentally measur- 
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able kinematical conditions, the method is useless. This is likely to be the case for 
the breakup of oxygen on nuclear targets, in order to extract the much needed in¬ 
formation on the relevant for the helium burning in massive stars pO 


For a comprehensive review on the relevance of the radiative capture reactions in 


astrophysics see the paper by Fowler [15 


The case of Coulomb breakup of as a tool to obtain the Be{p, j^B cross 
section at low energies has been investigated by many authors recently |^, |^, |^, 

In particular, in ref. ^0 it was shown that the nuclear contribution to the breakup of 
^B is negligible at forward angles below a threshold dictated by the onset of nuclear 
absorption and diffraction effects. This is a very useful result, since (in contrast to 
the breakup of oxygen) the much simpler semiclassical description of the Coulomb 
breakup maybe used. 

In this article we consider the Coulomb breakup of ^B at different bombarding 
energies, namely 5 MeV/nucl, 50 MeV/nucleon and 250 MeV/nucleon, e.g., at the 
Notre Dame University (USA), RIKEN (Japan) and the GSI (Germany). Radioactive 
beams of boron can be obtained at these energies at several other labs around the 
world. The advantage of performing a series of experiments at different bombarding 
energies is that the multipolarity content of the breakup can be disentangled with 
more accuracy. This is shown in section 2. 

The reaction ^Re(p, 'j^B is dominated by the electric dipole (El) multipolarity 
at very low energies, below the I"*" resonance at 633 keV. Higher order effects, or 
Coulomb reacceleration effects, act so as to distort the hnal relative energy of the 
fragments after the breakup. In ref. this problem was studied by means of a 
classical prescription, based on a Monte-Carlo calculation. Typel and Baur have 
recently studied the same problem within second order time-dependent theory. More 
recently, Esbensen et ah 0 have shown that while the effect is relevant for the 
breakup of projectiles on heavy targets, it is small for the breakup of ^^Be and 
large incident energies. 

It has been shown [IT, ^ that a correct treatment of reacceleration effects can 
only be accomplished by solving the time-dependent Schroedinger equation non- 
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perturbatively. This is because, the reacceleration is an effect related not only to 
the magnitude of the Coulomb breakup probabilities, which can be small. It is also 
an effect of distortion (caused by the Coulomb held of the target) on the relative mo¬ 
tion between the fragments after the breakup. The absence of such distortions are of 
vital importance to validate the Coulomb dissociation method. We report in section 
3 a calculation based on the numerical solution of the three-dimensional Schroedinger 
equation in a potential model for p-|-^i?e. The model used is similar to the one used 
by Esbensen et al. [^| for the breakup of and projectiles. An additional 
difficulty here is the inclusion of the Coulomb barrier for the p —^ Be relative motion. 
The results of this non-perturbative model are compared with previous calculations. 
It is shown that the reacceleration effect is much smaller than the experimental un¬ 
certainties that one achieves in the Coulomb dissociation experiments for the breakup 
of boron. The relative population of the s, p, d and / waves in the continuum is ob¬ 
tained. We also deduce the angular correlations, useful to separate the multipolarity 
content of the breakup process. As a byproduct of our calculations several predictions 
are made which are useful for future experiments. 

In section 4 we show that although the luminosity of electron beams are much 
higher than the ones obtainable in heavy ion facilities, it is hard to use this advantage 
to extract the photodissociation cross sections needed for astrophysics from inelastic 
electron scattering. However, the total yields of target breakup by inelastic electron 
scattering are not small for ideal experiments. Our conclusions are presented in 
section 5. 


2 First-order breakup of 

The extraction of the radiative capture cross section Be{p, l)^B from Coulomb 
breakup experiments depends on the validity of the first-order perturbation theory 
fim . When this approximation is valid the Coulomb breakup cross sections of ^B into 


5 




( 1 ) 


p Be in a continuum state with energy E is given by 

1 


d?ac 8 


rJ‘^71 

»)SUE), 


dE dn 5 (E + Q)3 

where E^ = E -h Q is the excitation energy, Q = 137 MeV is the binding energy 
of the system, and St,x{E) is the astrophysical S-factor for the electric, n = E, or 
magnetic, vr = M, multipolarity of order A. Detailed balance was used to link the 
ground state of {J'^ = 2+) to the proton (J'" = I"*") and the ^Be {J'" = | ) in the 
hnal state. The functions n^rxlEx, D) are known as the virtual, or equivalent, photon 
numbers. In ref. it was shown that one can calculate these quantities for all 
bombarding energies in the semiclassical approach including the effects of retardation 
and Coulomb repulsion. This is important since both retardation and the Coulomb 
recoil of the classical trajectory influence the values of nj^xiEx, at bombarding 
energies around 50-200 MeV/nucleon. At lower energies retardation plays no role 
and the results of ref. reproduce the semiclassical Coulomb excitation theory 


developed in ref. |^. At high energies, the Coulomb recoil is not relevant and the 
analytical results for the equivalent photon numbers are also reproduced [P^ . 

The differential cross section dac{E)/dE can also be written in the form of eq. (1), 
but with the virtual photon numbers integrated over the scattering angle 6, from zero 
to a maximum value corresponding to a grazing collision. To display the multipolarity 
dependence of the equivalent photon numbers (angle integrated) we plot them in hgure 
1 for the El, E2 and Ml multipolarities, and as a function of the excitation energy 
Ex- A lead target is assumed. We observe that the Ml multipolarity is about a 
factor 100, 10 and 3 smaller than the El multipolarity for bombarding energies of 5, 
50 and 250 MeV/nucleon, respectively. This is basically a consequence of the fact 
that umi oc {v/cYnEii where v is the projectile velocity. The E2 multipolarity is 
about 2-3 orders of magnitude larger than the El. This is a reflection of the “tidal” 
component of an 1/r^ force, which is more noticeable for small Fourier components 
of the held, and small bombarding energies. 

The astrophysical S-factors for the Be{p, y)®!? reaction is dominated by an 
resonance at E=633 keV. At lower energies the reaction is dominated by the El 
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capture from s and d waves in the continuum. A model calculation by Kim et al. 
yields an E2 cross section which is a factor 500-1000 smaller than the El cross 


section, over the relevant energy domain. In view of the results presented in hgure 
1 we conclude that the E2 breakup maybe as important as the El breakup mode in 
collisions with the bombarding energies studied here. This fact has been exploited 
in ref. to suggest that the data reported by Motobayashi et al. would 

imply an Sei-factor which is about 30% smaller than what was extracted from the 
experimental analysis. In view of the potentiality of the Coulomb dissociation method 
for this particular problem it is important to look at this question more closely. 

To calculate da/dE we will use the results of Kim et al. for the astrophysical 
S-factors of the El, E2 and Ml capture modes. The Ml and E2 capture modes are 
influenced by the continuum resonance at 633 keV. This influence is much stronger in 
the Ml case, so that this capture mode is only relevant for the energy region around 
633 keV. In hgure 2 we plot the energy spectrum of the Coulomb breakup of 
projectiles on lead at several energies and for the El, E2 and Ml dissociation modes. 
Apart from a small resonant component in the E2 dissociation mode, the El and E2 
contributions to the breakup are dominated by a prominent bump around 300-500 
keV in the relative energy of the fragments. This peak is the result of the folding 
of a increasing function of E, the Sommerfeld penetration factor exp{—27rp}, with a 
decreasing function of E, the equivalent photon numbers times the {E -|- Q)~^ factor 
in eq. (1). It is also clearly seen that the E2 component of the Coulomb breakup is 
comparable to the El component at the region of astrophysical interest {E <1 MeV). 
However, this result is strongly dependent on the angular kinematical conditions in 
measurements of the breakup cross sections, and also on the model for the capture 
cross sections used as input. 

To show this more clearly we plot in hgure 3 the angular distribution of the 
c.m. of the fragments with a relative energy of 100 keV, and for the El and E2 
multipolarities. A peculiar feature of the angular distribution is that the El breakup 
mode is concentrated at small scattering angles, while the E2 breakup mode is spread 
over a large angular region. This is more apparent in hgure 4, for E = 500 keV. 
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Eiab [MeV/nuch] 

<^{6 < 6m) [mb] 


a{6 < 6s) [mb] 



El 

E2 

El 

E2 

5 

572 

392 

159 

73 

50 

437 

210 

172 

23 

250 

178 

66 

38 

0.95 


Table I - El and E2 components in the Coulomb breakup of for events below a 
maximum scattering angle of 6m = 30°, 3° and 0.3° or 6s = 75°, 9.4° and 2.5° and at 
bombarding energies corresponding to 5, 50 and 250 MeV/nucleon, respectively. 


This important feature allows one to select events within an angular region and to 
manipulate the relative contribution of the E2 and El components to the breakup. For 
example, by selecting events for which the scattering occurs below a certain maximum 
angular value the E2 contribution can be reduced considerably. The calculations of 
the cross sections presented in figures 2-4 were obtained by integrating the equivalent 
photon numbers in eq. ( 1 ) from 6 * = 0 ° up to 0 = 6 * 5 , where 6s is the grazing 
scattering angle, above which the strong interaction sets in. These are 75°, 9.4°, and 
2.5° for Eiab = 5, 50 and 250 MeV/nucleon, respectively. Let us now assume that in a 
given experiment one selects events which occur below the angles of 30°, 3° and 0.3° 
for Eiab = 5, 50 and 250 MeV/nucleon, respectively. Then the total cross section, 
integrated over angle and energy, varies considerably and is shown in table 1. We 
notice that the E2 contribution is much more affected by the angular limitation of 
the events than the El component. 


We now compare our results with the experiment of Motobayashi et al. for the 
angular distribution of the breakup of on lead at 46.5 MeV/nucleon. We use the 
detection efficiency for the data with Erei = 600 keV obtained from the authors . 


In the experiment the uncertainty of the relative energy of the fragments amounted 
to 100 keV. Thus, we used as input in the calculations the S-factors averaged over 
energy bins of 100 keV and the model of Kim et al. However, for the El 

capture mode we multiplied the results of Kim et al. by a factor so that we could 
reproduce the data from ref. |^. After this procedure, the effective S'£;i-factor that 






























we use in our calculations is close to 18 eV b. Our results are compared with the 
data of Motobayashi et ah in the hgure 5. The solid line is the breakup cross 
section multiplied by the detection efficiency assuming only the El breakup mode. 
The dashed line includes the E2 breakup mode. In the experiment of Motobayashi et 
ah the restriction on the angular distribution of the fragments considerably decreased 
the effects of the E2 contribution to the breakup. But, if the model of Kim et ah 


2411 is valid they might not be negligible, as suggested in ref. and concluded 


from hgure 5. However, one must be skeptical about any conclusions related to the 
E2 contribution in this reaction. The model of Kim et ah has been criticised 
by Barker pj. In several other models [ 


191 the E2 and Ml contribution to the 


'^Be{p, 7 )®S reaction is considerably smaller than the ones predicted by the model 
of Kim et ah |^. If we use the model of Typel and Baur for the E2 capture 
model of the 'y)^B reaction we get the dotted line in hgure 5, which shows 

almost no inhuence of this multipolarity to the data of ref. |^. A recent discussion 
on the contribution of the E2-breakup to the data presented in hg. 5 has also been 


published in refs. |25 


We notice that the pattern of the theoretical angular distribution shown in hg. 5 is 
slightly diherent than the one shown in hg. 2 of ref. |^, particularly at small angles. 
The reason for these diherences is not quite understood, since we use the efficiency 


curve suplied by those authors. A possibility is that in ref. [T^ the efficiency curves 
were obtained by using a theoretical prediction for the dissociation modes as input to 
a Monte-Carlo simulation in which the detection efficiency and the geometry of the 
detectors were accounted for. This contrasts with the approach used here in which 
a simple folding of the theoretical prediction and the “model-dependent” efficiency 
curve is done. It has also been claimed ||^ that the angular averaging is not so simple 
as treated here, because of the poor experimental resolution and rather large energy 
bin sizes. It seems that this averaging, with proper inclusion of the detector geometry. 


has been made more accurately in ref. with the Monte-Carlo simulation. 

By switching off the retardation effects in the calculation of the virtual photon 
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numberg]^ which enter eq. (1) we obtain an average of 3 percent decrease of the cal¬ 
culated curves presented in figure 5, a small effect for this bombarding energy. But, 
these effects will increase with the bombarding energies and for higher excitation en¬ 
ergies, and should be considered in a more detailed analysis of future experiments. 
For the other (higher) relative energies presented in ref. |T^ we find that the ratio 
a{El + E2)/a{El) reduces considerably. In view of the discrepancies between the 
different nuclear models for the Be{p, lY'B reaction we feel that no quantitative 
conclusions can be drawn about the amount of the Ml and the E2 contribution to 
the Coulomb breakup of ^B. In fact, the Coulomb dissociation experiments will cer¬ 
tainly be very useful to solve this problem, since the several multipolarities contribute 
differently at different kinematical conditions. 

In figure 6 it is shown the energy spectrum of the breakup yields and compared 
with the experimental data of Motobayashi et ah |^. For qualitative comparisons 
we use the model of Kim et ah |Q. Again, we scale down the El capture cross section 
from Kim et al. [Q which amounts to use a constant factor Sei = 18 eV b. The 
solid curve in fig. 6 is our theoretical calculation (El only) folded with the detection 
efficiency. The dashed curve includes the E2 contribution to the breakup. In this 
case the same detection efficiency is used to obtain the E2 component. However, 
the E2 cross section was integrated only over the angular region where the detection 
efficiency (figure 2 of ref. |^) is relevant. As stated before (see figures 4 and 5), the 
exclusion of large angles, reduces considerably the E2 contribution to the breakup 
data presented in ref. [|T^]. Thus, it is not correct to assume the same detection 
efficiency for the E2 and the El breakup mode. If we assume that the El and E2 
efficiencies are the same, fig. 6 shows that only the point at 600 keV is appreciably 
modified by the inclusion of the E2 and of the Ml components. This is basically due 
to the presence of the I’*' resonance at 633 keV. 

We notice that the Sei =18 eV b value used in our calculations contrasts with 
that used {Sei = 15 eV b) in the calculation of ref. and presented in their fig. 

^This can be done by neglecting retardation effects in the calculation of the electromagnetic 
propagator, and on the calculation of the Coulomb trajectory M. 
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2. The reason for this difference might be due to their different numerical approach, 
as we discussed above in connection with the angular distribution presented in £g. 
5. However, we note that the S'd = 18 eV b adopted here seems consistent with the 
(very preliminary) value of Sei = 16.7±3.2 eV b given by those authors, and claimed 
to be the best £t to their full set of data. 

In hgure 7 we show the aE 2 /o'Ei ratio as a function of the scattering angle, the 
laboratory energy, and for relative energies between the fragments equal to 100 keV 
(a) and 500 keV (b). Once again, the model of Kim et ah m was used. One 
sees that the ratio decreases with the relative energy and with scattering angle. The 
ratio is also smaller for bombarding energies in the range of 10-200 MeV/nucleon. 
Experiments at different laboratory energies would allow a separation of the relative 
contributions of the El and E2 multipolarities to the breakup. 

As stated before, the validity of the Coulomb dissociation method to extract 
the astrophysical S-factors is strongly dependent on the validity of the hrst-order 
perturbation theory. In the next section we present a study of a non-perturbative 
approach to reacceleration effects. 


3 Non-perturbative Breakup of 

A non-perturbative approach to the study of reacceleration effects was proposed in 
refs. In the system of reference of we can write the time-dependent 

wavefunction for the p Be relative motion as 

T(r, t) = - Y, w;m(r, t) Yim{r) (2) 

r j 

Lm 

Initially the ^B is assumed to be in its bound state. We use a static p -1-'^ Be Woods- 
Saxon nuclear potential which reproduces the binding energy of ^B. The parameters 
used are Vq = —32.65 MeV, R = 2.95 fm and a = 0.52 fm for the depth, mean radius 
and diffusivity of the potential, respectively. The Coulomb potential of a uniform 
distribution of the '^Be charge with radius Rc = R is added. 
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We will treat the dynamical problem of the evolution of the wavefunction of the 
system with a time-dependent Coulomb interaction between the projectile and the 
target. However, we keep the nuclear potential constant. It is well known that a single 
set of potential parameters cannot reproduce simultaneously the binding, as well as 
continuum resonances of a system (see, e.g., the work of Tombrello |^). But, these 
are small changes in the potential parameters, and since we are interested here in the 
reacceleration caused by the time-dependent Coulomb interaction, we will neglect the 
dynamical modihcation of the nuclear potential. 

For simplicity we will assume that the projectile moves along a straight-line. To 
account for retardation effects we use the Lienard-Wiechert interaction between the 
target and the projectile |^. Expanding the interaction in lowest order at the region 
of the projectile we hnd that the time-dependent El potential is given by 



VEi{t) = \ — Ze'y r \ Fi _i(r) - Fii(r) b + \/2 v t Fio(r) \ f{t) , 


where 


1 - 3/2 


( 3 ) 


( 4 ) 


f{t) = p + 

In the equations above 7 = (1 — Z is the target charge, and b is the 

impact parameter. Eq. (3) is the effective El interaction. It is obtained by using the 
form (p0 — j.A/c^) for the electromagnetic interaction and the continuity equation 
to link the magnetic and the electric interactions. This is important for collisions 
where v/c is not negligible. The effects of retardation is manifest in the appearance 
of the factors 7 in eq. (3). We will restrict ourselves here to the El interaction. But 
our conclusions will be of general validity and will not be affected by the use of this 
approximation. 

When the expansion (2) and the nuclear plus the time-dependent interaction Vei 
is inserted in the time-dependent Schrodinger equation a set of coupled differential 
equations is obtained. The expansion is truncated at / = 4, for practical purposes, 
and the coupled equations is solved numerically. The numerical procedure is similar 


to that explained in ref. [17|. For example, if we restrict only to the time-dependent 
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El interaction, the equation to be solved is 


X 

+ 


/(/ + 1) 2/ifec jr / N 

- r '^<’'1 

I'm' 


fta(r. t)-^ m 

\i\( ‘ ^‘')-( ‘ 1 Ml 


\/2 v t 


I 1 V 
—m 0 m' 



t) 


“^f^bc duim 

h dt 


( 5 ) 


This equation couples uim with neighbouring angular states (due to the discrete values 
of the Wigner symbols) with A/ = ±1 and Am = ±1,0. Including the E2 interaction, 
complicates a bit more the coupled equations, but they are treatable with a truncation 
in the maximum /-value. 

We use a mesh of 1000 points for the coordinate r, with a mesh size of Ar = 0.1 fm. 
The time is discretized in a mesh of 800 points with a mesh size of 1 fm/c, starting 
at t = —200 fm/c. As explained in refs. 0,0, this numerical approach allows 
one to calculate the time dependence of the wavefunctions and observables useful 
to study the effects of reacceleration, e.g., occupation probabilities, momentum and 
energy shifts, etc. We will concentrate here on the results of this numerical approach 
which are relevant to the breakup. We will Erst assume a bombarding energy of 
50 MeV/nucleon. 

The time-dependent wavefunctions computed at 600 fm/c are not influenced by 
the Coulomb held of the target and are running waves (time-dependent continuum 
states) added to the remaining part of the ground state wavefunction (which is an 
I = 1 state). The continuum states are obtained by removing the ground-state part 
of the time-dependent wavefunction , i.e.. 




T(t)- < T(t) To > To 



< To|T(t) > 



( 6 ) 


In figure 8 the probability density of the p+’^ Be relative motion in the continuum 
states is plotted as a function of the radial distance, at t=400 fm/c, and for a collision 
with a lead target with an impact parameter 6=15 fm. Also shown is the ground state 
wavefunction. One sees that at this time the continuum wavefunction is already far 
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from the range of the static nuclear potential. It contains the characteristic behaviour 
of a spreading wavepacket. 

The wavefunction presented in hgure 8 is an admixture of several continuum an¬ 
gular momentum states. In hgure 9 we show the occupation probabilities of different 
continuum angular momentum states, for a collision with an impact parameter of 
6 = 15 fm and 6 = 50 fm, respectively. We observe that at 6 = 15 fm the states 
are more evenly occupied than at b=50 fm. This occurs because at b=15 fm the 
interaction is stronger and the continuum-continuum coupling is also stronger. The 
most relevant states are / = 0,1, 2, 3. At 6 = 50 fm the interaction is weaker, and the 
continuum-continuum coupling is small. In this case, hrst order perturbation theory 
should work well. In fact, at large impact parameters the most occupied states are 
the s-wave (/ = 0) and the d-wave (/ = 2) states. These are the same states which 


are relevant for the dipole capture from the continuum of the p +‘ Be system [24 


As stated above, of interest here is the comparison between the numerical results 
with the hrst-order perturbation theory. Only then we can access the magnitude of 
the reacceleration effect. In hrst-order perturbation theory the excitation amplitude 
is given by 




1 /■» 
ih J-c 


</14i(t) 


i > 


(7) 


where fiuj = {Ef — Ei) is the excitation energy. The time integrals can be done 
analytically. The amplitude is proportional to the electric dipole and quadrupole 
matrix-elements for a transition between the ground state and a continuum state 
with energy E = fiuj — Q. The continuum states are determined numerically for the 
same nuclear potential. We select the Z = 0 — 3 continuum states as input to calculate 
the matrix elements. 

To compare with the non-perturbative calculations, we Fourier transform the time- 
dependent wavefunction. One gets 

“ ( 8 ) 


where 


'J'c(p) = ^Cim{p) Yim{p) 

Im 


Cim{p) = Y “ J dr r ji{pr) Uim{r, t) 


(9) 
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where ji is the spherical Bessel function and uim{r, t) are the radial continuum wave- 
functions, with the ground-state subtracted. The probability density for an excitation 
to a hnal state with energy E is directly obtained from the above result. It is 


V{h, E, n) 


1 /2mbcN3/2 

2^ ) 


Ve 




.PJ 


( 10 ) 


Integrating over hi we get 


V{b, E) 


1 

2 


/2^x3/2 


Ve 


Vm(,P) 

Im 


( 11 ) 


The form of the energy spectrum of the relative motion of the fragments in hrst- 
order perturbation theory (solid) and in the non-perturbative approach (dashed) is 
shown in figure 10 . At small impact parameters the coupling between the continuum 
states is stronger and the reacceleration effect changes the form of the spectrum. Due 
to the reacceleration effect higher energy states are populated than in the perturbative 
calculation. At large impact parameters this effect is smaller, as clearly seen in £g. 


10 (b). 

The average energy of the relative motion can be compared with the results of the 
hrst-order perturbation theory. This can be obtained by using eq. ( 8 ) to obtain the 
first-order energy spectrum and the non-perturbative one described above. Averaging 
over impact parameter, the result is shown in figure 11 for the Pb reaction as 

a function of the bombarding energy (dashed line). 

In ref. a simple semiclassical model was used to describe the reacceleration 
effect. Neglecting the binding energy of the {a = b + c) system, the reacceleration 


causes a shift in the energy of the fragments which is given in the model of ref. 
by 


AEh = -AEr = 


71 


ZgZe^ 


4 


m, 




( 12 ) 


These energies are much smaller than the beam energy Eiab- The energy shifts can 
thus be calculated as small corrections to the velocities of the fragments. The modi- 
hcation of the hnal relative motion energy of the fragments due to the reacceleration 
effect is, in this approximation, equal to 


A 1 ^2 1 {AEbf 

2 4 Hbc 


( 13 ) 
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The results obtained by using this approximation is also shown in figure 11 (solid-line) 
after an average over the impact parameter (this average is weighted by the breakup 
probabilities, as explained in ref. |^). We observe that, despite the very different 
approaches used, the agreement between the two calculations is quite reasonable. 
Also shown in figure 11 is the Monte-Carlo method used in ref. to calculate the 
reacceleration energy from a classical prescription (solid circles). The error bars are 
inherent to numerical accuracy. We see that the numerical approach presented here 
yields the smallest result for the reacceleration energy (dashed-line). This can be 
understood as follows. In the approaches of refs. M the reacceleration energies 


were calculated by assuming that the breakup occurs at the distance of closest ap¬ 
proach. This overestimates the reacceleration energy, since in a quantum mechanical 
approach the proton is first brought to the continuum, then tunnels the Coulomb 
barrier, inducing a time-delay to the reacceleration process. This leads to an effective 
breakup position which is farther away than the distance of closest approach. This is 


more critical in the approach used in ref. [^. In any case, the reacceleration effect 
poses no major restriction to the Coulomb dissociation of ^B. 

In ref. it was shown that the reacceleration effects are very important for 
the breakup of projectiles on lead at an incident energy of 28 MeV/nucleon. 
However, it was also shown that the effect is negligible for the breakup of on 
lead at 72 MeV/nucleon. This result was ascribed to the higher binding of ^^Be and 
to the smaller effective dipole charge |^. For the breakup of ^B projectiles, the 
presence of the Coulomb barrier diminishes further the effect of reacceleration. Our 
results were obtained using a potential model for p Be, similar to that of Kim et 
ah 1^. As, discussed earlier, this model is thought to be very poor and unreliable 
1^. The use of a more refined model for the purpose of calculation of the effects of 
reacceleration is very difficult and beyond the scope of this work. 

Finally, we calculate the angular correlation of the fragments in the frame of 
reference of the p Be. This could be a useful tool to disentangle the El from the 
E2 dissociation mode in future experiments. In our model we can study this effect 
by varying the impact parameter, since as we have seen in figure 9, the occupation 
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probabilities for L=1 and 3 states (accessible also via direct E2 transitions) are also 
populated via the El interaction at small impact parameters. The angular correlation 
is just proportional to |T(p)p, obtained from eqs. (8) and (9). In hgure 12 we show 
the angular correlation of the fragments, in the frame of reference of the projectile, 
for E = 500 keV, 6^ = 4°, a bombarding energy of 50 MeV/nucleon, and 6 = 15 fm 
(dashed) and 6 = 50 fm (solid). Although small, the occupation of the 1=1, and 3 
states cause sizeable modihcation of the angular correlation. The simpler approach 
of Baur and Weber can be used due to the validity of the hrst-order perturbation 
theory. 


4 Photodissociation by Electron Beams 


The luminosity of electron beams are so much higher than the ones at heavy ion 
beam facilities that electron scattering experiments should in principle be another 
possible tool of obtaining the photodissociation cross sections useful for astrophysical 
purposes. While such procedure would be useless for the case, it could be useful 
for other radiative capture reactions of interest, e.g., the ^^C(a, 7)^®0 reaction. 

It is well known that for small momentum transfers electron scattering probes the 
same matrix elements as in photonuclear experiments (see, e.g, |Q). The condition 
is that the four-momentum transfer is much smaller than electron four-momentum 
(we use now h = c = 1), i.e., = {k' — k)^ k^. Also, for forward scattering, the 

longitudinal momentum transfer is 


{k' — k)s = k cos6* — — 


E' -E 


U! 


(14) 


where k^ is the component of the electron momentum along the beam axis. Thus, for 
the scattering at very forward angles, q = |k' — k| uj/c, i.e., the momentum transfer 
is the same as the one by a real photon. Under these circumstances, and in the long- 
wavelength approximation, both the transverse and the Coulomb matrix elements for 
inelastic electron scattering are proportional to the electromagnetic matrix elements 
It is easy to show that the differential cross section for electron scattering 
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becomes in this limit 


(15) 


d (Te _ 1 dflj^x T^X/ \ 

dVt doj ui ^ dVt 

where is the photonuclear cross section, and 

duMx a r£'E'- (Q.k)(Q.k') 


dfl 

duEx 

dVt 


TT^ 
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EE' + k.k' + m2 1 vEE' - (Q.k)(Q.kO - 
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where = 2EE'/2), and uj = E - E'. 

In order to guess the reaction yields due to the electron scattering at the “photon- 
point” {q = uj), we integrate the above equations over angles. Since they are a rapidly 
decreasing function of the scattering angle, we extend the integration to all angles 
to obtain an upper limit of the photon-point physics in electron accelerators. In this 
approximation the differential cross section dae/du can be expressed in the same form 
as eq. (15), but with the total virtual photon numbers given by 


^MA (l^) 

nEx{i^) 


2a 

7T 

2a 

TT 



A 1 

aTTJ ’ 


(17) 


where ^ = 2'yE/uj. 

There are two important differences between this result and the virtual photon 
numbers that one obtains from Coulomb excitation experiments. Firstly, the virtual 
photon numbers carry a factor due to the coherence effect of all charges in the 
nucleus. Secondly, the Compton wavelength of a nucleus is much smaller than that 
of an electron. The “smearing” of the electron wavefunction in the direction perpen¬ 
dicular to the beam weakens the strength of its interaction with the nucleus (see a 
discussion of this effect in ref. [^). Thus, the cross sections for photo-dissociation in 
electron scattering experiments are much smaller than in Coulomb excitation. This 
disadvantage could be compensated by means of a large beam luminosity. 

The cross section for the reaction ^®0(7, aY‘^C is about 3 nb for a relative energy 
of 1 MeV. Using the eqs. (18) we get ~ ~ 1 nb for CEBAF energies (U ~ 4 
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GeV). For a beam luminosity of 3 x 10^® cm?/s a considerable number of events can be 
obtained with a fixed target experiment. However, the extraction of very slow reaction 
products from the target poses a serious experimental difficulty. Targets of are 
impossible to construct due to the small lifetime (~ 700 ms). However, an attempt 
to use this method to study the would be worthwhile. As shown in ref. 

||34|| , coincidence experiments with the identihcation of the hnal products, i.e., the a 
and/or the oxygen, are quite possible. For astrophysical purposes a triple coincidence 
measurement is necessary. 

5 Conclusions 

The Coulomb dissociation of is a promising tool to extract the radiative capture 
cross section Be{p, 'j^B at low energies. The El and E2 contribution to the breakup 
are separable by using proper kinematical conditions. Instead of posing a difficulty, 
the enhancement of the E2 contribution to the Coulomb breakup at low energies 
is another useful tool to access information on the relative importance of the many 
multipolarities which might contribute to a certain radiative capture reaction, and 
which otherwise can only be extracted from dubious theoretical models. 

The effects of reacceleration are small and can be neglected. We have shown this 
with a series of different approaches, comparing perturbative and non-perturbative 
methods. 

I indebted to Drs. H. Esbensen, J. Bahcall and K. Langanke for useful discussions. 
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Figure Captions 


Fig. 1 The multipolarity dependence of the equivalent photon numbers (angle 
integrated) for the El (solid), E2 (dashed) and Ml (dotted) multipolarities, and as 
a function of the excitation energy E^. A lead target was assumed. Figure (a), (b) 
and (c) is for projectiles with bombarding energies of 5, 50 and 250 MeV/nucleon, 
respectively. 

Fig. 2 The energy spectrum of the Coulomb breakup of projectiles on lead 
at 5 (a), 50 (b) and 250 MeV/nucleon (c) and for the El, E2 and Ml dissociation 
modes. 


Fig. 3 Angular distribution of the c.m. of p Be with a relative energy of 100 
keV, and for the El and E2 multipolarities. Figure (a), (b) and (c) are for bombarding 
energies of 5, 50 and 250 MeV/nucleon, respectively. 


Fig. 4 Same as in fig. 3, but for E = 500 keV. 


Fig. 5 Angular distribution of the breakup of ^B on lead at 46.5 MeV/nucleon. 

m 


The data points are from Motobayashi et ah 


The solid line is a theoretical 


calculation for the El component folded with the detection efficiency. The dashed 
and dotted curves include the E2 component with different inputs of the E2 capture 
in the '^Be{p, 'y)^B reaction. 

Fig. 6 The energy spectrum of the breakup of ^B on lead at 46.5 MeV/nucleon. 

[0 


The solid line is a theoretical 


The data points are from Motobayashi et ah 
calculation for the El component folded with the detection efficiency. The dashed 
curve include the E2 component using the model of Kim et ah m as input for 
capture cross section '^Be{p, 'y)^B. 


Fig. 7 The ratio (Je 2 /o'ei as a function of the scattering angle and of the 
laboratory energy, and for relative energies between the fragments equal to 100 keV 
(a) and 500 keV (b). 
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Fig. 8 The probability density of the p Be relative motion in the continuum 
states as a function of the radial distance and at t=400 fm/c. The reaction + Pb 
at 50 MeV/nucleon is considered. 

Fig. 9 The occupation probabilities of s (1=0), p (1=1), d(l=2), f (1=3) and g 
(1=4) waves in the continuum, for a ^B Pb collision at 50 MeV/nucleon and an 
impact parameter of 6 = 15 fm (a) and 6 = 50 fm (b). 

Fig. 10 Energy spectrum of the relative motion of Be in first-order pertur¬ 
bation theory (solid line) and in the non-perturbative approach (dashed) at (a) b=15 
fm and (b) b=50 fm. A ^B Pb collision at 50 MeV/nucleon is assumed. 

Fig. 11 Total energy associated with the reacceleration effect in the Coulomb 
breakup ^B on lead targets, as a function of the laboratory energy per nucleon. 

Fig. 12 Angular correlation between the proton and the ^i?e-nucleus in their 
c.m. frame of reference in the Coulomb breakup of ^B on lead at 50 MeV/nucleon 
and a scattering angle of 4°. 
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